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Abstract 

The coincidence site lattices (CSLs) of prominent 4-dimensional lattices are considered. CSLs in 3 
dimensions have been used for decades to describe grain boundaries in crystals. Quasicrystals suggest 
to also look at CSLs in dimensions d > 3. Here, we discuss the CSLs of the root lattice A4 and the 
hypercubic lattices, which are of particular interest both from the mathematical and the crystallographic 
viewpoint. Quaternion algebras are used to derive their coincidence rotations and the CSLs. We make 
use of the fact that the CSLs can be linked to certain ideals and compute their indices, their multiplicities 
and encapsulate all this in generating functions in terms of Dirichlet series. In addition, we sketch how 
these results can be generalised for 4-dimensional Z-modules by discussing the icosian ring. 

1 Introduction 

Let us begin with a brief recapitulation of concepts and methods. 

1.1 Coincidence site lattices and modules 

Let r C M. d be a <i-dimensional lattice and Re 0(d) an isometry. Then, R is called a coincidence isometry 
of r if r(R) := r n Rr is a lattice of finite index in r, and r(R) is called an (ordinary or simple) 
coincidence site lattice (CSL) (for an introduction, see pQ). The group of all coincidence isometries of r is 
called OC (T) , whereas the subgroup of all orientation preserving isometries is called the group of coincidence 
rotations and referred to as SOC(r). The coincidence index S(R) is defined as the index of r(R) in r, and 
a = {E(R) I R £ OC(r)} denotes the (ordinary or simple) coincidence spectrum. These definitions can be 
generalised to include the possibility of multiple CSLs, see [21 OH]. 

All these concepts also extend to Z-modules [TJ[S]. By a Z-module of dimension d and rank k, we mean 
the Z-span of k rationally independent vectors in M. d . 

1.2 Quaternions and rotations 

Rotations in R 4 can be parameterised by a pair of quaternions [51 [7J [B] . Recall that the quaternion algebra 
H(M) is the vector space R 4 equipped with a special (associative but non-commutative) product. In standard 
notation, where l,i,j,k form an orthonormal basis of R 4 , the product satisfies the defining relational i 2 — 
j 2 = k 2 = — 1 and ijk = —kji = — 1. This definition is extended to all quaternions by linearity. All non-zero 
quaternions q = (a, b, c, d) = a + bi + cj + dk have an inverse q^ 1 = j^tQ, where q — (a, —6, — c, —d) is the 

conjugate of q and \q\ 2 = a 2 + b 2 + c 2 + d 2 is its (reduced) norm. With this notation at hand, we can write 
any rotation in R 4 as 

Rx = R(q,p)x = -r^—qxp. (1) 

\m 

** Corresponding author. Email: pzeincr@math.uni-bielcfcld.dc 

^^Note that these relations are very similar to the relations satisfied by the Pauli matrices. In fact, dividing the Pauli matrices 
by the imaginary unit gives a representation of A; by 2 X 2-matrices. 



Reflections and, more generally, all orientation reversing isometries may be parameterised by a pair of 
quaternions as well, but a conjugation is required in addition. For our purposes, it is sufficient to focus the 
discussion to (proper) rotations. This is no restriction, since all the lattices and modules we discuss have a 
reflection as a symmetry operation and thus every orientation reversing coincidence isometry can be written 
as the product of this particular reflection and a coincidence rotation with the same index. In other words, 
every orientation reversing coincidence isometry is symmetry related to a coincidence rotation with the same 
index (see below). Hence all CSLs are already obtained by (proper) coincidence rotations and the spectrum 
does not change either. There are exactly as many coincidence rotations as orientation reversing isometries, 
and the total number of coincidence isometries is just twice the number of coincidence rotations. 

1.3 Symmetry related rotations and CSLs 

Two coincidence rotations R and R' are called symmetry related if there is an element Q in the point group 
P(r) such that R' — RQ. Note that Rr = RT holds if and only if R and R' are symmetry related. 
In particular, two symmetry related rotations R and R' generate the same CSL, so that r(R) = r(R'). 
However, the converse statement is not true in general, i.e., two rotations generating the same CSL need not 
be symmetry related. In particular, in all the examples presented below, there are additional rotations that 
generate the same CSL. 

2 Hypercubic lattices 

The hypercubic lattices have been discussed in detail in fTJEJ. Here, we just review the most important facts, 
reformulate some of them and add some new details. In particular, we present an explicit expression for the 
number of CSLs. 

2.1 Centred hypercubic lattice 

The centred hypercubic lattice D|, which is the dual lattice of the root lattice D4, which in turn is a similar 
sublattice of D\ , can be identified with the Hurwitz ring J of integer quaternions. This is the Z-span (meaning 
the set of all integral linear combinations) of the quaternions (1, 0, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0), |(1, 1, 1, 1). Any 
coincidence rotation can be parameterised by a pair of integer quaternions. More precisely, a rotation is a 
coincidence rotation of D\ if and only if it can be parameterised by an admissible pair of quaternions [TJ [S] . 
Here, we call a quaternion q € J primitive if — q € I with n £ N implies n — 1. A primitive pair (q,p) is called 
admissible if \qp\ G Z. In addition, we call a primitive quaternion q € JJ reduced if \q\ 2 is odd (for a more 
general discussion of reduced quaternions, see [10] )■ Due to the unique (left or right) prime decomposition 
in J, we can decompose any primitive quaternion q £ J as q r s, where q r is reduced and |s| 2 is a power of 2. 
Note that q r is unique up to right multiplication by a unil|3 of JJ. 
For any admissible pair (q,p), we can define the factors 



The pair (q a ,Pa) ■— (& q q r , cipPr) is called the reduced extension of the pair (q,p). Note that two admissible 
pairs are symmetry related if and only if their reduced extension pairs are equal (up to units). The coincidence 
index of a rotation R(q,p) can now be written [TJ [5] as 



Since any (positive) integer can be written as a sum of four squares, this implies that the spectrum a is the 
set of all odd natural numbers, a — {1, 3, 5,7,.. .} 

The CSLs can be calculated explicitly [TTJ. They read 




and 




S Di {R{q,p)) = lcm(|g r | 2 , \p r \ 2 ) = \q a \ 2 = \p a \ 2 . 



(2) 




JJ 



+ 



Ip, 



'a- 



(3) 



2 Recall that the units are those quaternions u € J for which \u\ 2 = 1. 
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Since the point group of D4 contains 576 rotations, the number of coincidence rotations of a given index 
n can be written as 576 /j^ (n). This gives a multiplicative arithmetic function, i.e. 

f r n°t M = r D °l (m) fig (n) if m, n are coprime, (4) 

which is completely determined by /£, *(1) = 1, /£,°*(2 r ) = for r > 1 and 

/D 4 *b r ) = ^-yP r ~ 1 (p r+1 + P r_1 - 2 ) if P is an odd P rimc > r > !• ( 5 ) 
The multiplicativity of fjj*(n) suggests the use of a Dirichlet series as a generating function for it: 

D<u ,tt ^ s iia-^-xi-p 2 -) 

16 36 64 168 144 196 576 324 
= 1 + ^ + ^ + ^ + ^ + lP + lF + l^ + 17^ + '--- (6) 

Similarly, one can calculate the number f D (n) of different CSLs of a given index n. Clearly, one has f D (n) < 
f™*(n), and we do not have equality in general, since two CLSs that are not symmetry related may generate 
the same CSL. In fact, we have the following theorem which tells us when two CSLs are equal [TT] . 

Theorem 1 Let (qi,pi) and ((727P2) be two primitive reduced admissible pairs. Then, 

jn^l^jn^i (7 ) 
mpi\ mP2\ 

holds if and only if \q\pi\ = \q%P2\, lcm(|qi| 2 , |pi| 2 ) = lcm(|q 2 | 2 , \P2?), glcd(gi, \piq\\) = glcd(g 2 , |p2<?a|) and 
glcd(pi, \pm\) = glcd(p 2 , \p2q2D hold. 

Here, glcd denotes the greatest left common divisor in J, which is defined up to units of J. We find that 
frj 4 ( n ) is multiplicative, too, and is determined by /q (1) = 1, /£> 4 (2 r ) = for r > 1 and 

J^S (P 2r+1 + P 2r - 2 ~ , if r > 1 is odd, 

f Di (P ) \^ {p ^ +p 2r-2_ 2pr/ 2-^ if r >2iseven, 

for odd primes p. The corresponding Dirichlet series reads 

M _ V = TT 1 + P~" + 2 P'~ S + 2p- 2s + p 1 - 2 * + p 1 - 3 * 

dAS) ^ n s 11 (l-p 2 - s )(l-y- 2s ) 

n=l p^2 v f J\ / 

_ 16 36 64 152 144 196 576 324 
+ ¥ + 5 1 + 7 1 + ~W + IP + 13 7 + 15 7 + 17 7 + ' ' ' ' 

Differences to Eq. ([6]) occur for all integers that are divisible by the square of an odd prime. 
2.2 Primitive hypercubic lattice 

It is well known that the primitive hypercubic lattice Z 4 , which is a sublattice of D\ of index 2, has a smaller 
point group than D\, containing only 192 rotations, so that \Aut(D\) : Aut(l^)\ = 3. As a consequence, every 
class of symmetry related rotations splits into three classes, one of which has the same coincidence index 
as before, E^(R) = Sd 4 (R), while the other two classes have index 2£r) 4 (R). The number of coincidence 
rotations is thus given by 192/£J*(n), where /I°*(n) is again multiplicative, but slightly more complicated 
than (see [3[H])- It has the generating function 

$™* (s) = y &± = (1 + 2 1 - s )<F z ? t ( S ) = (1 + 2 1 -) TT ^ +P ~ S \ { n + P T\ 
n=l 

2 16 36 32 64 168 72 144 196 128 576 324 
+ 2 I+ 3 7+ 5 7+ 6 7+ 7 7 + _ 9 r+ HP + Tl 7+ 13 7 4 " 14 7 + 15 7 + T7 7 + " ' ' ^ 



3 



The CSLs split into two classes only, one of which has odd and the other even index. Hence, the generating 
function for the number of CSLs reads 



<Ms) = (1 + 2- s )$ D4 ( S ) = (1 + 2" s ) \\ 



l+p~ s + 2p'- s + 2p- 2s + p 1 - 23 + p 1 - 33 



(l-p 2 -)(l-pi- 2 *) 

( 1 16 36 16 64 152 36 144 196 64 576 324 
+ ~2 S + 3 1 + ¥ + 6 1 + 7 1 + ~W + IcF + IF + T3 7 + 14 7 + T5 7 + VP + ' ' ' ' 

Differences to Eq. © occur for all even integers and all integers that are divisible by the square of an odd 
prime. The spectrum is a = {2 e m I 6 {0, 1}, and m £ N odd}. 



3 Root lattice 

Usually, the A4 lattice — we will denote it by L in the following — is embedded in K 5 as a lattice plane. 
However, this is inconvenient for our purposes and we prefer to look at it in R 4 , since we want to exploit the 
useful parameterisation by quaternions, which we do not have at hand in 5 dimensions. A possible basis for 
L consists of the 4 vectors 

(l,0,0,0),i(-l,l,l,l),(0,-l,0,0),i(0,l,r-l,-r), (10) 

where r = 1+ 2 V ^ is the golden mean whose algebraic conjugate r 1 can be written as t' = — - = 1 — t. 
The lattice L cannot be identified with a ring of quaternions. However, if we interpret the basis vectors as 
quaternions, they relate to the icosian ring I, which is the Z[r]-span of the 4 quaternions 

(1, 0,0, 0), (0, 1,0, 0), ^(1, 1, 1, 1), ^(1 - r,r,0, 1). (11) 

Note that neither L nor I are invariant under algebraic conjugation. Combining the algebraic conjugation 
with a permutation of the last two components yields an involution of the second kind x := (x' , x' 1} x' 3 , x' 2 ), 
which was called twist map in [12[ 113] . Note that L — L is invariant under the twist map, which in addition 
is an antiautomorphism of I. The twist map is the key to our analysis since it provides us with a useful 
representation of L, 

L = {xel\x = x} = {x + x\xel} (12) 

and its CSLs — see below. 

Any coincidence rotation of L (and only those) can be parameterised by a single primitive admissible 
quaternion q G I by means of R(q) — j^qxq (see [IS])- Here, admissible means that \qq\ € N. 

Admissibility guarantees that one can define the factors 



~« ■ V s cd (M 2 > l«l a ) and a§ ' aq V g cd d«l 2 ' l«l 2 ) ' (13) 

which are elements of Z[r]. Note, however, that they are only defined up to units of Z[r]. In addition, we 
define the extension of a primitive admissible quaternion by q a = a q q. The coincidence index now reads [13] 

^ (*(?)) = A nf m = \l9\a q a s = lcm(| g | 2 , |q| 2 ). (14) 
gcd(|g| 2 , |g| 2 ) 

Consequently, the coincidence spectrum is N. An explicit expression for the CSLs exists, too: 

L n qaL ^ a = {q a x + xq a I x e 1} = (q a I + Iq a ) n L. (15) 

\QaQa I 

Along the same lines as for the hypercubic lattices, one can calculate the number of coincidence rotations 
and the number of CSLs of a given index n. The expressions obtained are slightly more complicated, since 
one has to distinguish between the three cases p — 5, p = ±1 (mod 5) and p = ±2 (mod 5), as these 
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primes behave differently and correspond to the ramifying, splitting and inert primes of Z[r], respectively. If 
120/^°*(n) is the number of coincidence rotations — recall that the point group of L contains 120 elements 
— its generating Dirichlet series is given by 

rt M = l+^ -pr (l+p-'Xl+p 1 -') pr l+pr* 
p=±l(5) P l\ Pi p = ±2 (5) F 

_ 5 10 20 30 50 50 80 90 150 144 
+ 2* + 3 1 + I 7 + 5 1 + 6 1 + 7 1 + ¥ + ¥ + W + IP + ' ' ' ' 

Similarly, we have for the number of CSLs 

5~ S \ -p-r 1 + p- S -,-r 1 + p- S + 2P 1 - 3 + 2 P - 2S + p 1 - 23 + p 1 - 33 



$a 4 (s) =1+6 



1 + 


p~ s 


n 


1 - 


p 2-s 


= ±1(5) 


50 


50 


80 


6 s 


+ 7"" 


h 8^ + 



7 p=±2(5) ■■- ■ " ' 



_ ! 5 10 20 6 50 50 80 90 30 144 

+ 2 T + ^ + 4 T+ 5 7 + ^ + F + ^ + ^ + l(F + Tl 7H ' 

Here, a criterion analogous to Theorem [1] exists, which will be published elsewhere [llj . 



4 Icosian ring 

It is also interesting to discuss the coincidences of the icosian ring I itself. Note that it is no lattice in 4-space 
but a Z-module of rank 8. However, this does not matter and one can argue similarly to the hypercubic case. 
The main difference is that one has to work with Z[r] instead of Z, which does not cause any problems. In 
fact, some steps are even easier since the prime 2 does not play a special role here. So, a primitive quaternion 
is already reduced and hence the notion of reducibility is not needed here. 

Let us call a pair (q,p) £ I x I primitive admissible if q,p are primitive and \qp\ S Z[r]. Then, a rotation 
is a coincidence rotation if and only if it is parameterised by a primitive admissible pair {q,p). For any 
primitive admissible pair (q,p), we can once more define 



^ and ~ - / ^ 



y g cd(M 2 , |p| 2 ) y gcd(M 2 ,w 2 r 

which are again elements of Z[r], defined up to units. As before, (q a ,p a ) := (a q q,a p p) is the corresponding 
extension pair. We get the following expression for the coincidence index |llj : 

El(R(q,p)) = N(lcm(|g| 2 , |p| 2 )) = N(| fe | 2 ) = N(|p Q | 2 ). 

Here N(a) = \aa'\ is the (number theoretic) norm of a € Z[r]. Thus, the coincidence spectrum consists of all 
integers that contain prime factors p = ±2 (mod 5) only with even power. The CSMs read explicitly 

flip 

inf-^ = q a I + Ip a . 

\1P\ 

If 7200/]f o *(n) denotes the number of coincidence rotations and fi(n) the number of CSMs, we get again 
f{ ot (n) > fi(n). Again, equality fails to hold in general, and an analogue of Theorem Q] tells us which CSMs 
are equal [TT] . The corresponding Dirichlet series read 

t (i+s-xi+s 1 -) -,-r / (l+p-'xi+p^n 2 n (i+p- 2s )(i+ P 2 - 2s ) 

+ 4 7+ 5 7 + _ 9 r + TF + lo 7 "^ 19 7+ 20 7 + 25 T + '" ' ' 
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and 

n s )~ (l-5 2 - s )(l-5 1 " 2s ) 11 b (1 -p 2 - ;s )(l -p 1 - 28 ) 

1 + p~ 2s + 2p 2 - 2s + 2p- 4s + p 2 - 4s + p 2 - 6s 

X 



n 



(l-p4-2«)(l _p2-4«) 
p=±2(5) V 1 n F ' 

25 36 100 288 410 400 900 912 

= H 1 1 1 1 1 1 1 1 • (17) 

4 s 5 s 9 s 11 s 16 s 19 s 20 s 25 s v ' 
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